
■ S37 





THE 



ONTARIO WATER RESOURCES 



COMMISSION 




SIMULATED THERMAL EFFLUENT 



LAKE ONTARIO 



APRIL, 1968 



Copyright Provisions and Restrictions on Copying: 

This Ontario Ministry of the Environment work is protected by Crown copyright 
(unless otherwise indicated), which is held by the Queen's Printer for Ontario. It 
may be reproduced for non-commercial purposes if credit is given and Crown 
copyright is acknowledged. 

It may not be reproduced, in all or in part, for any commercial purpose except 
under a licence from the Queen's Printer for Ontario. 

For information on reproducing Government of Ontario works, please contact 
ServiceOntario Publications at copyright(5>ontario.ca 



SIMULATED THERMAL EFFLUENT 
MTO LAKE ONTARIO 



APRIL, 1968 



M. D, Palmer, 

Water Quality Surveys Branch,, 

Division of Sanitary Engineering 



WQTS-1 






SIMULATED THERMAL EFFLUENT 
INTO LAKE ONTARIO ' 

ABSTRACT 

The dispersion pattern of Duffin Creek water into Lake 

Ontario was determined under two different wind conditions in the 
spring time when the flow was 18 cfs and 5 F warmer than the lake 
water. Continuous dye injection was used to trace the plume. Based 
upon experimental results,,, a model was developed which predicts maximum 

concentrations within 25 percent for distances up to 40! 00 feet from source. 

— 8 
Maximum dye concentrations varied as the (distance) ' * . The model 

applied to an existing thermal generating station at Lakeview was 

capable of predicting maximum, temperatures within 25 percent of 

measured values. 



SIMULATED THERMAL EFFLUENT 
INTO LAKE ONTARIO 

INTRODUCTION 

There is relatively little information available that allows one 
to predict what effect an effluent released at the shore will have on 
adjacent areas. Obviously, this type of information is mandatory for the 
proper locating of municipal service intakes and outlets as well as the 
effect of industry on adjacent areas, The information must be formulated 
in such, a way that it can foe used for different loading and longshore current 
conditions. The Duffin Creek study was designed with the equations 
formulated having three variables; loadings, longshore currents and 
distance from the source. In addition, the equations incorporate a 
temperature difference between the effluent and the receiving body which 
is likely to be the case for most discharge conditions. The equations 
were developed from controlled experiments in the field, then checked by 
applying them to a completely different area. They were also compared 
to results from similar studies on Lake Huron as a test of their general 
application. While the equations developed are obviously not universal 
m capable of handling all physical situations, they do provide a valuable 
prediction vehicle for small discharges at the shore. As more field studies 
are conducted, it will be possible to expand the equations to incorporate 
more physical factors and produce more accurate predictions. 
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An Interim approach to the problem would be to check the 
application of the prediction equations with limited field surveys , then 
decide whether it is necessary to conduct full studies or accept the ex- 
isting equations* In each case, the results should be such that they 
can be extended and maximized by utilizing longterm meteorological data, 
OUTLINE OF PROBLEM 

At pre sent , a thermal generating plant is under construction at 
Pickering. The plant on completion in 1971, will discharge warm cooling 
water (10-25°F above lake temperature) at a rate of 1200 cfs. A study was 
conducted at the mouth of Duff in Creek to arrive at an estimate of the 
dispersion pattern. This stream which is tributary to Lake Ontario was 
selected because of its geographical proximity to the plant, the lack of 
a well formed delta, and the large marsh just upstream of the mouth which 

ponds and routes the flow. To simulate a thermal effluent, the study 

o 
was conducted in the spring-time when the flow was 180 cfs and 5 F 

warmer than the lake water. Before conducting the study, a mathematical 

model was developed for the analysis of the dispersion data. 

MATHEMATICAL MODEL 

One form of the steady state one dimensional dispersion 

equation is: 

£ ' 6 T& W 
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where C = dye concentration ppm 

U = mean velocity fps 

^C = distance ft. 
The constant " £ " called the dlffusivity has been found to be a variable. 
Diffusivity ( ^ ) Is a function of mean velocity gradients, turbulence and 
thermal stratification with the result that various techniques have been 
designed to define the variability of vT , Some of the most successful 
methods in natural stream flows have utilized the measurement of mean 
velocities in two dimensions (Fischer/ 1967), However, these methods are 
not easily adapted to dispersion in the lakes where velocities along the 

shore are 0.1 **~G,3 fps. Csanady in the Bale du Dore study 

(Csanady, 1964) found diffusivity to be highly variable in both dimensions 
in Lake Huron, 

G, Gifford (1955, 1957, 1960) proposed a dispersion model 
for a continuous source based upon a 2- dimensional Gaussian distribution 
using measured values of the variances in two dimensions. This 
method uses a mass balance criteria. J.E. Foxworthy et al (1966) 
successfully applied Gifford 1 s model to the dispersion study at Los Angeles, 
for a point source . 
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mass flow of dye in lbs/sec. 

* 2 

average variances it . 
depth measured from the surface in ft. 
distance measured parallel to the water 
surface perpendicular to SA ft. 
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Equations (2) and (3) are for a point source and require 
modification for a steady state area source as follows (See Foxworthy 
(1966) for detailed development and assumptions): 
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The variances were determined on an individual point basis 
as suggested by Foxworthy 
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A short bibliography on lake currents is included to assist in 
the selection "U" values for computations without drogue studies for 
estimating purposes. 
EXPERIMENT 

Rhodamine B dye was injected (see Fig. 1) at a rate of approxi- 
mately 200 ppb (a flow gauge existed two miles upstream of outlet). The 
distribution of the dye was not uniform at the outlet due to the non- uniform 
mean velocity (see Fig. 4), The surface layer and the eastern side of 
the creek were actually flowing upstream due to the wind conditions. 
Consequently, the mass flow term " Q" was determined by integrating the 
concentration and velocity profiles at the outlet. The longshore currents 
were measured by tracking drogues at a two foot depth placed 5 00 feet 
directly south of the outlet while wind velocities were measured with an 
anemometer, (Table 1). Once the direction of the current was determined, 
marker buoys were established at the various range lines (constant dis- 
tances from creek outlet) using an optical range finder. Dye was 
continuously injected at a steady rate. Once the plume had reached 
approximately 6000 feet from the source, sampling was started on the 
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nearest range line. The boat traversed the dye plume sampling at a 
depth of two feet to determine the location of the maximum concentration 
in the "Z" direction. At this point, samples were taken at one foot 
intervals (y-direction) to determine the depth to the maximum concen- 
tration. The vessel then made a complete traverse from well off-shore 
to the shore, sampling at the depth of maximum concentration. 
Positions were determined from the shore using a transit mounted 
optical range finder with a one meter base. All concentrations were 
measured with a fluorometer. 

Dye was normally injected for about three hours with the 
result that the wind conditions varied little (see Table 1). Two runs 
were made under different wind conditions resulting in a Western and 
Eastern dispersion pattern, but in each case, the temperature of the 
creek water was 44°F. and the lake temperature 500 feet offshore was 39~F, 

RESULTS 

In both studies, the maximum concentration of dye remained 
very near to the shoreline at each range line. Sampling was restricted 
to depths of approximately 1.5 feet due to the draft of the small boat. 
It was observed that the dye plumes were wider dispersed at 500 feet 
than at 1000 feet. This is probably due to the inertia of the creek flow 
which is directly south with a mean velocity of 0.5 fps. The inertia of the 
creek water must be dissipated, then redirected along the shore to 
conform with the wind driven lake currents along the shore. Beyond 
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this initial spreading at the mouth, the concentrations on the average varied 
as x in the westerly direction and as x * in the easterly direction. 

Csanady (1967) found the concentration varied approximately between 
x~ 1#0 and x" 0,7 . The confinement of the thermal regime produced higher 
concentrations in a narrower zone along the shore. 

The dye concentration profiles at the various distances from 

— 2 
the mouth are presented in Figs, 2,3,4,5 and 6 and the variances ( Q" ) 

for each profile appear in Fig , 7 » The distribution of dye with depth is 
fairly uniform and appears to have a constant variance of approximately 
one foot (Fig, 6), This is probably due to the bottom having a uniform 
awa y- fro m- shore gradient , 

To apply the prediction model, it was necessary to develop a 
relationship between the variance and distance from the outlet, Foxworthy 
fitted a straight line to a plot of log O" vs log x. However, this method 
produces erroneous results when applied to Duff in Creek due to the 
initially large plume at 5 00 feet from the creek mouth, A two segment 
line was fitted to the measured variances by joining the variance at the 
creek outlet and the variance at 1000 feet, then a second line was fitted 
for distances greater than 1000 feet. This method produced a good pre- 
diction model for maximum concentration for Duffin Creek, but failed to 
define the initial large spreading of the plume at the creek outlet which 
would probably be dependant on the creek flow and cross section. This 
study also failed to define the relationships between variances and different 



thermal lakewise temperature gradients, which are likely a factor as 

— 1 — 7 

Gsanady' s concentrations vary approximately as x to x "in the 

-G • 8 
summer and the Duffin Creek study vary approximately as x 

Further, the study failed to define the relationship between variances 

and longshore currents. 

MODEL COMPUTATIONS - WESTERN PLUME 

The model maximum concentrations were computed by sub- 
stituting the following values into equations 5 and 6: 
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X 





U 


CTu 

ft* 


°» 


feet 


lbs/ sec 


fps 


ft 





11.1 x 10" 4 


0.58 


3.16 


5 .36 


500 


n.i x io" 4 


0.43 


0.94 


36 


1000 


.-4 
11.1 x 10 


0.33 


0.97 


42 


2000 


11,1 x 10" 4 


0.33 


0.97 


67 


4000 


11.1 x 1Q" 4 


0.33 


0.97 


108 



The variation in velocity downstream was accounted for by employing 
dimensionless velocity ratios to equation 5. The results are as 



follows: 
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Distance 


' C 

max 

Measured 


C 
max 

Model 


Percent 
Difference 


ft 


ppb 


ppb 


% 


WESTERN PLUME 








500 


130* 


150 


18 


1000 


205 


170 


17 


2000 


130 


112 


14 


4000 


72 


70 




EASTERN PLUME 








500 


140 


175 


25 


1000 


100 


94 


6 


1500 


60 


70 : 


13 


1 






; 



True maximum missed. 



A temperature and drogue study conducted at an existing 
thermal generating station* (Lakevlew ) with a cooling water flow of 
1160 cfs at a temperature of 62.6°F. In May showed that the warmer 
water stays near the surface and the shore (see Fig. 9). The sub- 
surface cooling water had acquired lake temperature by 4000 feet 

.- - o 

while the surface cooling water Is still 4 F. warmer than the lake water 

at 4000 feet. The Duff In Creek dye plumes were similar In configuration 
to the measured Lakevlew thermal plumes. 
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The Duff in Greek model was applied to the Lakeview thermal 
plant site to predict maximum downstream temperatures, based on 
temperature and flow data at the outlet (assuming constant along shore 
currents and heat as a complete conservative with no atmospheric heat 
transfer) the following results were obtained: 



Distance from 


Predicted Maximum 


Measured Maxi- 


Difference 


the Outlet 


Temperature 


mum Temperature 




feet 


o f 


°F 


/a 


500 


60.5 


62.6 


7 


1000 


57.0 


62.6 


18 


2 000 


53.4 


60.8 


25 


4000 


50.7 


53.4 


12 



It must be appreciated that the transposed model is an 
estimate which can only be applied to along shore plumes. A more 
comprehensive model would take into consideration current regimes, 
thermal stratification and atmospheric transfer and probably wave 
con siderations . 
CONCLUSIONS 

A mass flow model based on measured concentration pro- 
files for continuous constant dye injection was developed for a creek 

with a flow of 180 efs into Lake Ontario in the presence of a thermal 

o 
difference of 5 F, between the creek and the lake. The model predicts 
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dye concentrations within 25 percent of measured values in the effluent 

— ft 
plume for a mile from the outlet and vary as the (distance) u '°, The 

developed model was then used to predict the maximum temperatures in 
the effluent plume of an existing thermal generating station and was 
capable of predicting maximum, temperatures within 25 percent of the mea- 
sured values. 
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TABiE 1 



WIND CONDITIONS 

PUFFIN CREEK - APRIL 10/68 



Time 


Force (Knots) 


Direction 


Activity 


WESTERN 








0900 


4 


o 

245 




1000 


7 


280° 


dye injection 


1115 


10 


285° 




1330 


10 


280° 


sampling 


15 00 


11 


280° 




15 40 


7 


280° 




EASTERN 








07 00 


5 


110° 




0815 


5 


1 1 0° 




0930 


4 


95° 




1110 


m 


95° 




1150 


Q 

© 


113° 




1330 


15 


60° 




1430 


11 


5 0° 


dye injection 


1700 


1 


50° 




1815 





0° 


sampling 

i 



LIST OF SYMBOLS 
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APPENDIX 



x - distance measured in the direction of the dye 

plume movement - feet 
y - depth measured from the surface - feet 

z - distance measured parallel to the water surface 

and perpendicular fed x - -feet 
U - mean velocity in x direction fps 



~" -•* 2 

O variances in y and z direction ft 



- mass flow rate of dye lbs/sec 

fl - frequency 

c - concentration ppm 
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EASTERN PLUME 
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VERTICAL DYE CONCENTRATIONS 
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VARIANCE <*£ VS. DISTANCE 








LEGEND 
• WESTERN PLUME 








©EASTERN PLUME 






200 


M 


















100 


* 


• 






50 
i 


n 


° y 


o 

• 




J 30 

i 










H"20 










10 
B 

! 2 

1 

5 5 

> 


Ik 






i 







i S 8 


1 




10 


100 500 1000 


^000 


4000 




DISTANCE FROM OUTLET (FT) l» 




, 



fig. a 

MODEL AND EXPERIMENTAL 
RESULTS 



t 



250 
200- 

150- 



m 

0. 
ft. 

Z 

O 100 



< 
or 

i- 

2 

ilil 
O 

z 
o 
o 

itj 

> 
a 



50 



x 
2 



id 




LEGEND 

— -• WESTERN PLUME MEASURED VALUES 

— — ©EASTERN PLUME MEASURED VALUES 

▲ WESTERN PLUME THEORETICAL VALUES 

& EASTERN PLUME THEORETICAL VALUES 



X™ 



— X WESTERN PLUME THEORETICAL STRAIGHT 
LINE FIT 



1 






100 



500 1000 

I STANCE FROM OUTLET (FT) 



2:000 



4000 






*Tt E 13yDD,DDD c il73^ 



